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SUMMARY 

In an earlier analysis the influence of peripheral heat conduction was considered for 
the condition in which only the broad walls of rectangular ducts were conducting. The 
heat was generated solely within the broad walls, and the narrow side walls were non- 
conducting as well as unheated. Both the flow and the heat transfer were assumed to be 
fully developed. 

In this report the previous work has been extended to include two cases of greater 
generality. The first case was a rectangular duct of arbitrary aspect ratio with uniform 
heat generation and heat conduction within the broad walls as in the previous analysis. 

The extension made here is to allow the short walls to also be heat conducting, but with 
the restriction that a junction of infinite thermal resistance exists between the broad and 
short sides so that heat flow around the corners cannot occur within the walls. 

The temperatures within the duct were governed by a Poisson equation, with the nor- 
mal derivatives at the boundaries expressed in terms of the second partial derivatives 
along the boundary to account for the conduction within the walls. An analytical solution 
was obtained for the temperatures in the fluid and on the boundaries by using a product 
solution in the form of a Fourier series expansion. The unknown wall temperatures were 
contained in the integrals required for the Fourier coefficients so that the wall tempera- 
tures were expressed by an integral equation. An analytical solution was found, and re- 
sults were evaluated numerically on a digital computer. 

The peripheral conduction within the side walls produced the following effects: 

(1) The temperature distribution along the unheated side walls became quite uniform 
for even small values of wall conductivity. 

(2) For aspect ratios less than 5, the temperatures of the heated broad walls were 
lowered by a significant amount; for aspect ratios greater than 5, only small temperature 
decreases occurred on the broad walls for a very narrow region adjacent to the corner. 

The second case studied was a square duct with uniform heating and peripheral heat 
conduction in all walls. The solutions provided quantitative information on how the wall 
temperatures become more uniform as a function of wall conductivity. For large wall 
conductivity, the results were in agreement with the uniform peripheral wall tempera- 
ture solution. 


INTRODUCTION 


In reference 1 an analysis was made of the influence of peripheral wall conduction on 
the temperature distributions around the perimeter of rectangular channels that were 
cooled by fully developed laminar forced convection. The channels were heated along 
part or the entire width of only the broad sides with the short sides left unheated. Wall 
heat conduction was assumed to take place only within the broad sides; the short sides 
were of insulating material. Solutions were evaluated for aspect ratios from 1 to 00 with 
the wall conductivity varying from 0 to °°. In the present report the analysis of refer- 
ence 1 is extended to two more general cases in which all the channel walls are heat con- 
ducting. These cases are described as follows. 

Case I: The cases considered in reference 1 were models of typical cooling passages 
employed in flat-plate fuel assemblies of nuclear reactors. In each assembly a series of 
parallel heat generating plates with coolant flowing between them are supported by un- 
heated side plates. In general, all walls including the side plates are heat conducting 
with some heat flowing from the fuel plates around the corners into the side plates. Be- 
cause of the complexity of the most general problem, it was assumed in reference 1 that 
the side plates were nonconducting. In this report some cases of greater generality are 

treated that include the side wall conduc- 
tion. An attempt was made to account for 
the heat flow around the corner within the 
walls, but the results for that problem 
were unsatisfactory because of poor con- 
vergence of the Fourier series represent- 
ing the wall temperature derivatives at the 
corners. In the present analysis, there- 
fore, it was assumed that no heat flow oc- 
curred around the corners within the walls. 
This assumption eliminated some terms 
from the more general solution, and the 
convergence in this case was quite good. 
The situation with no heat flow around the 
corner represents a limiting case some- 
times encountered in practice where a 
poor joint exists between the fuel plates 
and the side support plates, thus providing 
a high thermal resistance at the corners. 
Figure 1(a) illustrates this case. The 
heat generation, which is either supplied 
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thermal 
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^ Peripheral conduction in both 
broad and short sides 

(a) Rectangular channel heated only on broad sides, with peripheral 
conduction in all sides and no heat flow around corner within the 
walls. 


Uniform heating 

id ± i 
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^Peripheral conduction 
in all sides 

(b) Square channel with uniform heating 
and peripheral conduction in all sides. 

Figure 1. - Models of channels analyzed. 
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to or generated within the walls, is taken to be 
uniformly supplied over the broad sides. The 
solution was obtained and the temperature distri- 
butions were evaluated for aspect ratios of 1, 2, 
5, and 10 for values of the wall conductivity pa- 
rameter that ranged from 0 to high values. 

Case II: The second case considered the 
effect of peripheral wall conduction on the tem- 
perature distribution in a square duct where all 
four walls conduct heat and are uniformly heated. 
This geometry is one that is sometimes con- 
sidered for honeycomb type fuel assemblies in advanced nuclear propulsion devices. The 
duct configuration is shown in figure 1(b). Temperature distributions were obtained for 
several wall conductivities between 0 and a high value. 

For both cases the fluid motion and the convective heat transfer are assumed to be 
laminar and fully developed with constant fluid properties. The same type of analytical 
approach as that in reference 1 is used. An additional part is incorporated into the so- 
lution to account for heat being conducted within the two short walls assumed to be insu- 
lated in the earlier paper. The additional part leads to a coupling of the temperature 
distributions for adjacent sides and a simultaneous solution is required. 

A similar analysis for turbulent convective heat transfer is not possible because the 
present state of knowledge of the turbulent mechanism is incomplete, especially for rec- 
tangular ducts. The analysis for laminar convection is of value because there are some 
conditions under which the laminar convection does occur, during reactor shutdown for 
example. Secondly, the solution for laminar convection will demonstrate in a qualitative 
way the influence of peripheral conduction on the wall temperatures. The laminar solu- 
tions would be expected to provide larger variations of wall temperature around the duct 
perimeter as compared with the turbulent case. 



Figure 2. - Coordinate system for model of rectangular thin- 
wall channel employed in analysis. 


ANALYSIS 


A coordinate system for rectangular (including square) ducts is defined as shown in 
figure 2. The momentum equation governing the velocity distribution under the condi- 
tions of laminar fully developed flow is given by 


*>= Ji? h + 

dz \3x 2 



( 1 ) 
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(Symbols are defined in appendix A. ) When heat is being transferred to the fluid, the 
temperatures throughout the fluid are determined by the energy equation which, under the 
assumed fully developed heat-transfer condition, takes the form 


= + + < 2 ) 
p V 2 ^ W 

where viscous dissipation is assumed negligible. Since the fluid properties p, c , k^, 
and ju are considered to be constant throughout the flow field, the momentum equation (1) 
and the energy equation (2) are not coupled, and the velocity profile can be determined 
from equation (1) directly. This has been done and the profile available in Knudsen and 
Katz (ref. 2, p. 101) is (after rearrangement into a more convenient form) 



2 _ 128 y 1 /l - e- n7Ty > 

3 v _5 Z-f n 5L _-n Try 

yn n=l, 3,5... n \ 1 + e j 


The energy equation can be simplified by making a heat balance for a volume extending 
over the entire channel cross section and having an incremental length dz. This heat 
balance results in the following relation: 


- 0T b 

4abpc u = Q 

p Sz 

where Q is the total heat input per unit channel length and is a constant for all the 
cases considered here. For fully developed conditions the temperature profiles have the 
same shape for all z so that 


3T _ 0T b 
0z 3z 


Hence, 
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Equation (4), Poisson’s equation, is to be solved using the velocity profile given by equa- 
tion (3) and is subject to the particular boundary conditions of interest in this analysis. 



y - b 


Figure 3. - Terms for heat balances at boundaries with peripheral conduction in walls (for unit 
axial length of duct). 
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Case I - Rectangular ducts with uniform internal heating throughout the broad walls 
only, and with peripheral heat conduction in all walls (no heat conduction around the cor- 
ner within the walls) . - The boundary conditions for this class of problems are derived 
by performing heat balances on elements of the heated and unheated walls. The terms 
entering into the heat balances are shown in figure 3 (p. 5). Note that only one quadrant 
of the duct cross section need be considered because of the symmetry that exists about 
the coordinate axes. The wall thickness w is considered to be sufficiently thin so that 
the temperature is constant through the wall thickness and is equal to the local fluid 
temperature at the wall. The local internal heat generation per unit area q^ will be 
held constant along the width of the heated wall, although the analytical technique used 
here can be applied for any prescribed variation of q^ with y. The total heat input per 

unit channel length in the flow direction Q is related to q by the equation Q = 4bq . 

& & 

Under these restrictions the thermal boundary conditions for the duct quadrant are as 
follows: 


x=a, 0<y<b 


0<x<a, y = b 


x = 0, 0<y<b 


0<x<a, y = 0 


de _ 1 , w jSv a 2 e 
ax b kf ay 2 

de = w k w s 2 e 
3 y" kf ax 2 


ie = ol 

3x 


M = o 
3y J 


Symmetry 

conditions 


(5a) 


(5b) 


(5c) 

(5d) 


The second derivative terms in (5a) and (5b) account for the peripheral conduction 
within the walls. When the wall conductivity k^ is 0, the problem reduces to one of 
uniform heat flux transferred to the fluid from the broad walls of rectangular channels 
with the short walls nonconducting, a set of conditions that has been treated in refer- 
ences 3 and 4. 

One other condition must be specified to complete the mathematical description of 
the problem. The condition is imposed that no heat be conducted around the corner within 
the walls at x = a, y = b, and this is expressed by equating the derivatives along each 
wall to 0 at the corner. Within the short wall, y = b, this condition gives 
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at x = a 


(6a) 


0X 

In the broad wall, x = a, the condition is 

— = 0 at y = b (6b) 

3y 

An analytical solution will now be obtained for the energy equation (4) subject to boundary 
conditions (5) and (6). 

The solution is derived by applying the method of superposition. As discussed in 
reference 5 (p. 9), if for a given geometry a solution of Poisson's equation is known for 
simple boundary conditions, then more complicated boundary conditions can be accounted 
for by superposing solutions of Laplace's equation. For the present problem the unknown 
dimensionless temperature 0 is replaced by a sum of a Poisson solution and two Laplace 
solutions: 


e - 0 pI + 0 L1 + e h2 


( 7 ) 


Equation (7) is substituted into the energy equation (4) and the boundary conditions (5) 
and (6). The resulting relations are divided as follows into three distinct partial differ- 
ential equations and their boundary conditions: 


v V. = ± ifei! 

ab u 


(8a) 


0<x<a, y = 0, b 


00 


PI 


3y 


= o 


(8b) 


x=0, 0<y<b 


00 


PI 


0X 


= 0 


(8c) 


x = a, 0<y<b 


00 


PI 1 


0x b 


V 0 L i = ° 


(8d) 

(9a) 
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00 L1 

0<x<a, y - 0 JjA - 0 

dy 

(9b) 

0<x<a, y = b ^ = w 

** *t Sx 2 

(9c) 

x-0, a, 0<y<b ljl -0 

3x 

(9d) 

y20 L2 = 0 

(10a) 

dO-r 2 

0 < x < a, y = 0, b = 0 

9y 

(10b) 

30j o 

x-0, 0<y<b ^ - 0 

9x 

(10c) 

n . , 90 L2 \/ 0 2 0 

x-a, 0<y<b - w 

9X k f 0y 2 

(10d) 


Poisson’s equation (8a) and boundary conditions (8b), (8c), and (8d) account for the 
total heat addition to the channel. This is necessary because a solution to a Laplace 
problem cannot add or subtract any net amount of heat to or from a closed region. The 
solution for 0^1 adjusts the fluid temperature distribution to account for a redistribu- 
tion of heat along the short side wall due to the presence of conduction within the wall. 
The 0^2 solution accounts for the conduction along the width of the broad wall. Thus, 
the original problem has been divided into a forced convection problem with simple 
boundary conditions and two problems of steady- state conduction. 

The solution for is derived in reference 4, and for convenience the final re- 
” A 2 2 

suits are summarized in appendix B. Both v = 0 and V Q-^2 ~ ® can sa tisfied 
by product solutions of the same general form: 


'Ll 


OO 

z 


B n cos 


/n7rx 


n=l, 2, 3* 


Va 


Hff) 


(ID 
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and 


oo 


'L2 


z 


j=l, 2, 3- 


Aj cos cosh 


( 12 ) 


where equations (11) and (12) each satisfy their zero normal derivative boundary condi- 
tions (9b), (9d), (10b), and (10c). To satisfy the remaining boundary conditions, (9c) 
and (lOd), equations (11) and (12) are each differentiated once (eq. (11) with respect to y 
and eq. (12) with respect to x), substituted into their respective boundary conditions 
(9c) and (lOd), and then expanded in a Fourier series to give 




(13) 



(14) 


When each integral on the right side of equations (13) and (14) is integrated twice by 
parts, and conditions (6a) and (6b) applied in the process, the Fourier coefficients Aj 
and B n are expressed as 


and 



(15) 


(16) 
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The coefficients in equations (15) and (16) are then substituted into equations (11) and 
(12), and the resulting expressions are placed in dimensionless form. The hyperbolic 
functions are converted to exponential form for greater convenience in numerical evalu- 
ation, and 0 L1 and 0 L g become 


'Ll 


oo 

I 


n=l, 2, 3- 


- 2nffK cos(n7rX) 


,-n 7 ry(l-Y) + e -n7iy(l+Y) 


1 - e 


-2n7ry 


/ 0(X, l)cos(n7rX)dx| 


(17) 


OO 



3-1) 2, 3- • • 


-H(l-X) -Jl(l+X) 


2]ttK 


COS(j7TY) 


+ e 


_ 2jn 
1 - e y 


[jT 1 e(l,Y)cos(j7rY)dYj 


(18) 


The solution for 0 is obtained by combining equations (17) and (18) with equation (7): 


oo 

- £ 


e(X,Y) = 0 PI (X,Y) - ) 2mrK cos(n?rX) 

n=l, 2, 3 


e -n77y(l-Y) + e -n7ry(l+Y) 


1 - e 


■2n7jy 


f 1 0(X, l)cos(n7rX)dX 


C0S(jffY) 

y 


3=1, 2,3- 


-1?(1-X) -il(l+X) 


+ e 


2 j 7T 

1 - e y 


Jjf 1 e(l,Y)cos(j7rY)dYj 


(19) 


If equation (19) is evaluated at X = 1 with Y variable and at Y = 1 with X variable, 
the result is two simultaneous integral equations for the unknown temperature distribu- 
tion 0(1, Y) along the heated wall and the temperature distribution 0(X, 1) along the un- 
heated wall: 


10 



oo 

-Z 


0(l,Y) = 0 pi (l,Y) - > 2mtK(- l) n ] 

n=l, 2, 3 


-mjy(l-Y) + e -n7ry(l+Y) 


1 - e 


-2mry 


K e(x ’ 


l)cos(n7TX)dXj 



^ 37r ^ cos(jnY) 

y 


J=l, 2,3- 


_ 2j7T 
1 + e y 


_ 2j77 

1 - e y 


EC e(1 - 


Y)cos(jnY)dYJ (20) 


oo 

-Z 


8(X, 1) = 0 pi (X, 1) - ) 2mK cos(nTTX) 

n=l, 2, 3 


1 + e ” 2niry 


1 - e 


-2niry 


f 1 6(X, l)cos(n7rX)dX 
Jc\ 



2]ttK. 

r 


j=l, 2, 3- 




(l-x) -21(i+x) 


+ e 


_ 2jj 

1 - e y 


f 1 d( l,Y)cos(jjrY)dY 

-O' 


( 21 ) 


These integral equations can be solved by a method given in reference 6. 

The integrals / 6(X, l)cos(n?rX)dX and f 0(1, Y)cos(jirY)dY are constants for 

*'0 T) 

each value of the summation indices n and j. These can then be defined as 


C = f 1 8(X, l)cos(n7rX)dX 

II 


( 22 ) 


D. s f 1 0( l,Y)cos(jirY)dY 
3 J 0 


(23) 


These constants are unknown as they contain the unknown wall temperatures 0(X, 1) and 
0(1, Y). Equation (20) is now multiplied by cos(mjfY)dY and integrated from Y=0 
to 1: 
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J 1 0(l,Y)cos(mirY)dY = f 1 fl pi (l, Y)cos(mirY)dY 


^ 2nrfK(-l) n C n J ^ n7?y(1 ' Y) + 
n=l, 2, 3* • * ^ 


e -n7iy(l+Y) 

2njry 


cos(m?rY)dY 


_ 2m7r> 

mffK / 1 + e ^ _ 

_2mrr * m 

1 - e y 


(24) 


In the evaluation of the last term on the right the following identity was employed: 


l 


cos(j7rY)cos(mjrY)dY = - when j = m 

0 2 


= 0 when j ^ m 


The term on the left side of equation (24) is equal to D . The first term on the 


right side is defined as 


E m 2 pi (l,Y)cos(mirY)dY 


(25) 


The integral in the second term on the right is carried out, and equation (24) can then be 
written as 


D = E 
m m 


OO 

I 

n=l, 2, 3- 


2m7i\ 


2n7rK(-l) n C - (.i. 1 )” 1 ™? 1 D (26) 

H n n 1 e\ ni 


{niryy + (m7r) 


2 r 


2m7r 


a - e 


Equation (26) is solved for D m and the result written in the form 
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where 


D e + 7 f C 

m = m Z-/ mn n 

n=l, 2,3- •• 


( 27 ) 



(28) 


m, lX n 


-2K(-1) (-1) 


mn 


7 



1 + (-) 2 


(29) 


In a similar fashion when equation (21) is multiplied by cos(lirX) and integrated the fol- 
lowing relation is obtained: 


C £ = S £ + £ h «j D j 

j=l, 2,3- - - 

where 




1 + £7tK 


1 + e 

ll-e 


-2liry 

-2Stny 


J 1 0 PI (X ’ 1 ) C0S (^ 7rX ) dx 


(30) 


(31) 


(32) 
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( 33 ) 


-2K(-l) £ (-l) j 

7 \lr 


4 + e “ 2£,7r \ 

l+f77K( 1 + e ] 

/ Gv 

1 + (— ) 

\l - e- 21 ’ 1 '/ 

\i / 


Equations (27) and (30) are the sought after simultaneous relations for the unknown inte- 
grals and D m that are contained in the solution for 0(X, Y) (eq. (19)). If equa- 
tions (27) and (30) are combined, a single expression is obtained for either the C^’s or 


00 oo 


D m = l e m + E Sn'mn) + E E < f mn h ni> D l < 34) 


n=l, 2, 3- • • / n=l, 2,3--- 1=1, 2,3- •• 


or 


oo oo 


= ( g l + E e j h £j) + Z ^ v * , lj 1 jn 

j=l, 2, 3* - * / j=l, 2, 3* • • n=l, 2, 3- 


E <M in >C n (35) 


Each of these relations is an infinite set of simultaneous equations for either the D m or 
C^. A truncated set of these equations was solved on a digital computer by the technique 
given by Kantorovich and Krylov (ref. 5, p. 20). Generally about 100 D m and 
were required to have the series solution for 6 converge. The and C„ decreased 

at least as well as l/rn and 1/ i . The quantities E m and that appear in the co- 
efficients e m and g^ were evaluated from the 0pj solution and their analytical forms 
are given in appendix C. The solution for 6 is now complete; however, to be usable in 
physical problems it must be given relative to the fluid bulk temperature. The dimen- 
sionless bulk temperature 0^ is defined as 
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The 0pj k is known from reference 4 and is given in appendix B. Rather than evaluate 
0 L1 b anci ^L2 b by direct integration, a much more simple approach is obtained by 
applying Green’s second identity. The procedure will be shown for only the function 
®L 1 k as the result for 0^ ^ is found in a similar fashion. The u/u in the b 

integral is eliminated by use of equation (8a), which gives 


0 


Ll, b 


=/ 7 ' 


e L1 v 2 e P1 


dx dy 


(37) 


The second form of Green's identity relates the volume and surface integrals for a 
volume. For the present problem, the velocity and temperature profiles are fully de- 
veloped and we can choose as a volume a unit length of the duct. Green's identity then 
provides the following relation composed of a boundary integral (not a line or contour 
integral) and an area integral over the cross section: 



(38) 


where d/dv denotes the normal derivative, ds the elemental length on the boundary, 

2 

and T the boundary. Because V = 0 (eq. (9a)), the second integral vanishes. The 
first integral contains normal derivatives of 0pj and 0pj on the region boundaries and 
these are 0 except at the channel walls. Hence equation (38) takes the form 

dX (39) 

1 


e Ll.b = ^ 


00 


pi 


sx 


x= 


f e L 1 (^)dY-i / 


®pi(X, 1) 


90 


Ll 


0Y 


Y= 


From equation (8d) 00pj/0X| = 1/y, and the 0j_^ an( j 30^/ 3Y| ^ ^ are evaluated 

from equation (17). These quantities are substituted into equation (39) and integrated to 


give the final result: 


0 Ll,b 



n=l, 2, 3- • - 


2KC. 


(njr) 2 F 


n 



(40) 


Similarly, 


15 


I 


(41) 


®L2, b = £ 7 <i,|2D i E i 

j=l, 2,3--- 

In summary, the final solution for the dimensionless temperature distribution 
throughout the fluid cross section is given by 

(o ~ 0 h ) = (0 pi - 0 pi b ) + (e L1 - 0 L1) b ) + (0 L2 - 0 l 2, b^ 

where (0pj - 0pj b ) is given in appendix B, 0^ by equation (17), 0^ by equation (18), 
0Li b by equation (40), 0^ b by equation (41), and the integrals D m and C7 in the 
Fourier coefficients are found from equations (34) and (35). 

Case II - Square duct with uniform internal heating and peripheral heat conduction 
throughout all the walls . - The analysis for this case will be presented in less detail than 
that for case I as the basic method is the same. If the dimension 2a is used for the 
length of the duct side, the relation between the heat generation per unit wall area and 
the total heat input per unit channel length is 


Q = 8aq g 

The boundary conditions analogous to those in equation (5) then become 


0 < y < a 

30 _ 

— + 


3 2 0 

(43 a) 


3x 

2a 

k t 

3y 2 


< a, y = a 

30 = 

:J- + 


3 2 0 

(43b) 


3y 

2a 

k t 

3x 2 


x = 0, 0 < 

y < a 

30 ; 

= 0 


(43 c) 



3x 




0 < x < a, 

o 

II 

30 

= 0 


(43d) 



3y 





From the symmetry of the problem there is no heat flow within the walls around the 
corner; hence, the derivatives at the corner remain 0 as in equations (6). 


i 
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The solution is again formed by a superposition of three solutions as in equation (7). 
The 0pjj (the use of II in the subscript distinguishes this solution from 0 pI in case I) 
solution must now account for heat addition along all sides of the duct so that the condi- 
tions in equations (8) now become 


2„ _ 1 u(x,y) 

v 0 pn ~ — — — 

a 2 u 


(44a) 


x = a, 0<y<a 


90 


PII 1 


9x 2a 


(44b) 


0<x<a, y = a 


90 


PII 1 


9y 2a 


(44 c) 


x = 0, 0<y<a 


90 


PII 


0X 


= 0 


(44d) 


0<x<a, y = 0 


90 


PII 


3y 


= o 


(44e) 


The conditions for 0^ and 6 ^ remain the same as in equations (9) and (10) except 
that the dimension b becomes a for the square configuration. 

The solution for 0 p jj is found from appendix B by letting y = 1 and /3 = 1, the 
latter value corresponding to the condition of uniform heating all around the duct 
boundary. 

The solutions for 0^j and 0^ can be obtained from equations (17) and (18) by 
letting y = 1: 


0 L1 = 


oo 

Z 


2n?rK cos(n7rX) 


r mr(l-Y) + e -mr(l+Y) 


1 - e 


-2n7T 


'L2 


n=l, 2, 3- 


OO 

Z 

3=1, 2,3- 


l( 


2)ttK cosQttY) 


r j7T(l-X) + e -jrr(l+X) 


1 - e 


•2jjr 


0(X, l)cos(n7rX)dX (45) 


] 


j^ 1 0(l,Y)cos(j7rY)dY] (46) 
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Since for a square duct these expressions are symmetric in X and Y, 


f 1 0(X, l)cos(n 7 TX)dX = f 1 0( l,Y)cos(n7TY)dY s S 

•'O J o 


( 47 ) 


When the S n are found, then 0^ and 0^ will be known and can be added to 0pjj to 
provide the temperature distribution. The S n are obtained in the same manner as the 
C n and Dj in case I. However, since the C n and would be equal for the present 
symmetric conditions, equations (27) and (30) simplify to 


S = p + 
m *m 


v S 
mn n 


n=l, 2, 3* 


(48) 


where 


m 


m 


1 + (m7rK) I 


1 + e 


>2m7r\ 


^1 - e 


■2mi t i 


(49) 


V 

mn 



(50) 


and 


P m = J e pn (X, l)cos(mwX)dX = 0pjj(l, Y)cos(m7TY)dY (51) 

Analytical expressions for P m are given in appendix C. Equation (48) is an infinite set 
of simultaneous equations for S m and was solved as in case I by the methods given in 
reference 5. Generally about 200 S m were required to obtain accurate wall tempera- 
ture distributions. 

To be applicable physically the solution has to be related to the bulk temperature 
0^. For the 0 piI function the bulk temperature is found from 0 p ^ in appendix B by 
letting y = 1 and (3 = 1. The values 0 L1 b and 9^ which are equal to each other, 
are found by following the same procedure leading to equation (40). The only change is 
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that 90pjj/3X | is now 1/2 rather than 1 fy as used for equation (40) and the no- 
tation changes as C n and F n are replaced by S n and P n , respectively. This yields 



In summary, the temperature distribution is given by 

6 ' 9 b = ( 6 PU " e PH, b) + pLl " fl Ll,b) + ( 0 L2 " 6 L2, b) ^ 53) 

The Poisson solutions 0 p jj and 0 p jj b are found from appendix B by letting y = 1 and 
/3 = 1. The Laplace solutions 0 L ^ and 0^ are given in equations (45) and (46) where 
the integrals equal to S n are found from the solution of the simultaneous equations (48). 
The bulk temperatures 0 ^ ^ and 0^ ^ are obtained from equation (52). 

RESULTS AND DISCUSSION 
Case I 

The preceding analysis was used to evaluate the dimensionless temperature distri- 
butions (1) throughout the fluid cross section for a square duct (figs. 4(a) and (b)) and 
(2) along the walls for ducts of various aspect ratios (figs. 5(a) to (d)). The isotherms 
of figures 4(a) and (b) were drawn through temperatures that were calculated from equa- 
tion (B2) with y = 1 and (3 = 0 for figure 4(a) and from equation (42) for figure 4(b). In 
figures 5(a) to (d), for the dimensionless wall temperatures, the origin of the abscissa 
corresponds to the center of the heated sides. As the curves are followed to the right, 
the temperature distributions proceed along the heated sides to the corners (from Y = 0 
to Y = 1) and then along the unheated sides from the corners to the center of the un- 
heated sides (from X = 1 to X = 0). The dimensionless temperatures are governed by 
two parameters - the channel aspect ratio y and the parameter K, which is referred to 
hereafter as the wall to fluid conductivity parameter. Each part of figure 5 is for a dif- 
ferent aspect ratio, y = 1, 2, 5, and 10. The solutions from the present analysis are 
given by solid lines, and for each solid curve the corresponding solution of reference 1, 
where only the heated walls are heat conducting, is given by the dotted lines. The dot- 
dash lines provide the limiting case when the walls are all nonconducting. 
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T ' T b 



Figure 4. - Isotherms and approximate heat flow directions within fluid for a square duct uniformly heated in two 
opposite sides only. 


An examination and comparison of the isotherms of figure 4, and of the dotted and 
solid sets of curves in figure 5, clearly reveal the influence of having the unheated walls 
become heat conducting. Two effects are evident. First, the temperatures along the 
unheated walls become much more uniform. Secondly, there is a reduction in the tem- 
peratures of the heated walls. This reduction is appreciable in channels of small aspect 
ratio (y = 1 and 2) and becomes smaller when y is increased to 5 or 10. To under- 
stand this behavior consider that in a square or rectangular duct the convection is poor 
in the corners because of the low fluid velocities. When only the heated walls are made 
conducting some of the heat generated near the corner will be conducted peripherally 
within the walls away from the corners to a region of better convection. This peripheral 
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T ' T b 

Q/4k f 



Figure 4. - Concluded. 

wall conduction thereby reduces the corner temperatures. Once the heat is redistributed 
over the conducting walls, it must all be transferred directly from these walls to the 
fluid. This is indicated in figure 4(a) which, although it is for nonconducting walls 
(K = 0), is typical of all square ducts with heating and conduction in opposite walls only. 
This is because conduction in only the two heated walls does not appreciably change the 
wall temperature distributions as evidenced by the dotted curves in figure 5(a). 

When the unheated side walls are also allowed to be heat conducting, these walls 
provide an additional path for the heat to flow from the heated walls into the fluid. This 
is seen by examining the patterns of the isothermal lines within the fluid and the resultant 
heat flow directions for the case when y - 1 and K = 5, as shown in figure 4(b). Some 
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Y * y/b X - x/a Y - y/b X ■ x/a 


(a) Aspect ratio, y * 1. (b) Aspect ratio, y - 2. 

Figure 5. - Wall temperatures with peripheral conduction in (1) no walls, (2) heated walls only, or (3) all walls with no heat 
flow around corners. 

of the heat generated in the heated walls now flows through the fluid in the corner regions 
to a portion of the unheated walls adjacent to the corners. This heat is then conducted 
along the unheated walls and is transferred to the fluid in a region of improved convec- 
tion. Hence, the unheated walls now become somewhat effective as heat-transfer sur- 
faces, and the temperatures along the heated walls are thereby reduced. 

One other feature of figure 4(b) is worth noting because it will aid the understanding 
of the wall temperature curves in figure 5. The heat transfer across the corners causes 
the isotherms to curve rapidly and group together in the corners. This produces very 
rapid temperature changes along the walls near the corners and obscures the fact in 
some cases that the derivatives along the walls are 0 at the corners as imposed on the 
analysis by equations (6). 
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(c) Aspect ratio, y - 5. (d) Aspect ratio, y = 10. 


Fiqure 5. - Conch irfpri 


As the aspect ratio is increased, the fact that the short wall occupies less of the 
total channel perimeter makes the conduction in the short wall proportionately less ef- 
fective in lowering the temperatures on the heated wall. For an aspect ratio of 10 the 
conduction in the short unheated walls produces only a small reduction in the tempera- 
tures of the heated walls close to the corner. 

The effect of increasing the parameter K is also demonstrated in figure 5. As K 
increases, the temperatures along all walls become more uniform, but have a higher 
value on the heated walls than on the unheated walls. It is this temperature difference 
that provides the driving potential for the heat flow from the heated walls, across the 
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fluid in the corner regions to the unheated 
walls. To obtain an indication of the mag- 
nitudes of the K parameter for some 
physical conditions, a channel is consid- 
ered where the thickness of the walls is 
equal to 1/20 the spacing between the 
broad heated walls, w/2a=l/20. For 
aluminum walls with water as a coolant, 

K is approximately 30. If the wall ma- 
terial is changed to stainless steel. K is 
reduced to about 3. If the channel is 
stainless steel and the coolant is a liquid 
metal, K becomes less than 1. 


Case II 


Figure 6. - Wall temperatures of square duct with uniform 
heating and peripheral conduction in all walls. 


The wall temperatures presented here 
are for a square duct that has uniform 
heating and peripheral conduction in all 
walls. Because of symmetry the tempera- 
ture distributions are the same for all 
sides; hence, only the distributions ex- 
tending from the center of one side to the 
corner need be given. Results are shown 

in figure 6 for K from 0 to 75. For a small K of unity there is a substantial reduction 
in the peak temperature and the temperature gradients along the wall. For a K of 25 
or greater the temperatures have become uniform and should correspond to the limiting 
case of a duct having a uniform peripheral temperature with uniform heat addition in the 
axial direction. As given in reference 3, the Nusselt number for this boundary condition 
in a square duct is 3. 60. When the definitions 


Nu = 


h2a 


and 


q 


g 


= _Q 

8a 


h < T w - T b> 
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are used, there is obtained 



Q/4k { Nu 3.60 


The value 0. 278 is in agreement with the limit for large K in figure 6 and also serves 
as a check on the correctness of the solutions. 

CONCLUSIONS 

5> ‘i 7 

A previous analysis on the effect of peripheral wall heat conduction on the convective 
heat transfer in rectangular ducts has been extended to study two situations of greater 
generality and complexity. The peripheral heat conduction within each wall was formu- 
lated in terms of an integral equation which was coupled with the energy equation in the 
fluid. The solutions were obtained in analytical form and evaluated for various aspect 
ratios and wall conductivities. 

Case I: A rectangular duct was considered with uniform heating and peripheral con- 
duction in the broad walls. When the short unheated walls were allowed to conduct heat 
but no heat was allowed to flow around the corners within the walls, the following changes 
were found as compared with a duct with nonconducting short walls: 

1. For moderate values of the wall conductivity the temperature distribution along 
the unheated wall became almost uniform. 

2. For duct aspect ratios less than 5 the conduction in the unheated wall lowered the 
dimensionless temperatures along the heated wall. The decrease in the difference be- 
tween dimensionless wall and bulk temperatures was from 15 to 30 percent depending on 
the conductivity and aspect ratio. For larger aspect ratios the temperatures along the 
heated wall were decreased only a small amount for a region near the corner. 

Case II: A square duct was considered, heat conducting and uniformly heated on all 
walls. Moderate values of wall conductivity were found to produce substantial reductions 
in the peak temperatures that occurred at the duct corners. For high values of the wall 
conductivity, the duct approached the known case of constant peripheral wall temperature 
with uniform heat addition in the flow direction. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 24, 1965. 
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APPENDIX A 


SYMBOLS 


A 

cross-sectional area 

k 

thermal conductivity 

5 i 

Fourier coefficients defined by 

Nu 

Nusselt number 

eq. (16) 

P m 

integral defined by eq. (51) 

A „ 

a 

Fourier coefficients defined by 
eq. (B7) 

half-length of short side of rec- 

P 

Pm 

static pressure 
coefficient defined by eq. (49) 


tangular duct or half-length of 
side of square duct 

Q 

total heat-transfer rate to fluid 
per unit channel length 

B n 

Fourier coefficients defined by 

q B 

heat flux through broa<^valls 

eq. (B9) 

q g 

heat generation rate in wall per 

B n 

Fourier coefficients defined by 

unit wall area 

eq. (15) 

q S 

heat flux through short walls 

b 

half-length of broad side of rec- 

S n 

integral defined by eq. (47) 


tangular duct 

s 

element of surface area around 

C n 

integral defined by eq. (22) 


duct perimeter per unit duct 

C P 

specific heat of fluid at constant 


length 


pressure 

T 

temperature 

D s 

integral defined by eq. (23) 

u 

local fluid velocity 

E 

m 

integral defined by eq. (25) 

u 

integrated mean fluid velocity 

e m 

coefficient defined by eq. (28) 

v mn 

coefficient defined by eq. (50) 

F f 

integral defined by eq. (32) 

w 

thickness of channel wall 

f mn 

coefficient defined by eq. (29) 

X 

dimensionless coordinate, x/a 

G 

coefficient defined by eq. (B4) 

X 

coordinate measured from center 


coefficient defined by eq. (31) 


of duct in direction parallel to 
short sides 

h 

heat-transfer coefficient 

Y 

dimensionless coordinate, y/b 

h £j 

coefficient defined by eq. (33) 

y 

coordinate measured from center 

K 

wall-to-fluid conduction parame- 
ter, wk^akj 


of duct in direction parallel to 
broad sides 


26 



z coordinate measured along 

channel length 

/3 ratio of heat fluxes, q^/q^ 

y channel aspect ratio, b/a 

9 dimensionless temperature, 

4k f T/Q 

p absolute fluid viscosity 

p fluid density 


Subscripts: 

b integrated bulk mean value 

f fluid 

L Laplace solution 

P Poisson solution 

w wall 

I case I 

II case II 

1, 2 distinguishes the two different 

Laplace solutions 
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APPENDIX B 


SOLUTION FOR 0 pi AND 8 PII 

The solution for Poisson’s equation 0 p - 0p ^ given in reference 4 is for rectangu- 
lar ducts that are being uniformly heated on each of the four walls, but where the heat 
flux qg on the short sides is an arbitrary fraction 0, between 0 and 1, of the flux q p 
on the broad walls. That is, 


q s = 0q B (Bl) 

When 0 = 0, all the heating is transferred from only the broad walls, and the ^ 

solution needed for case I results. When 0=1, the heat flux is a constant all around 
the duct periphery. If in addition to 0 = 1 the aspect ratio is set equal to 1, the 
e pn - 0pjj b solution for case II is found. 

The general Poisson solution is 

°P ‘ 0 P,b = e ir + e * + d l + e 2 " 0 P,b 


where each of the 6 terms on the right side is given as follows: 



m=l, 3, 5- • • 


0*(X, Y) = * (X 2 - y 2 Y 2 ) (B5) 

2 (y + 0) 
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e 2 ( x > Y ) 


z 


B cosh(n^Y) cos(njrX) 
sinh(nTry) 


(B6) 


(B7) 


(B8) 
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APPENDIX C 


EXPRESSIONS FOR E m , F f , AND P m 


oo 



3 = 1 , 2 , 3 


(j7iy) 2 + (n nr) 2 


(Cl) 



(C2) 

The coefficient P m is given by either of the following two forms resulting from the two 
integrals in equation (51). Although they both give the same numerical results, two dif- 
ferent algebraic expressions are found because 0 p jj has not been used in a symmetric 
form: 


P 


m 


-(-l) m 

(m7r) 2 


'G + _1\ eoth(m7r) 
s8 2/ 2 



I 

1=1 9 


(-l)j(-l) m jTT 
(jjr) 2 + (mu) 2 



(C3) 


or 
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r 


00 


= G j 

(-llMW-e] 4(-i> m \ 

1 

[ 1 + w tanh /pm! I 

i 

24(m7r)^ 7r® / . 4 

-*-gf 

l 2 ' 2 '\\ 


S. p=l, 3, 5- • • 

L \2/J 

J 


+ i^ + cotMm I I (Bm) y < C 4) 

2(m77)^ y ^ (rTr)^ + (mn-)^ y 


r=l,2,3- 


The Fourier coefficients A m , (A m ) , B^, and (B.) 
and (B9) and G is given by equation (B4). 


are given by equations (B7) 
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